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ABSTRACT 


Local and average Nusselt numbers have been experimentally 
determined around the circumference of a uniformly heated 
Pa eciteeiMuaemcylinder placed near a plane surface in a 
cross flow of air. Experiments were carried out for Reynolds 
mambers ranging trom 16,000 to 153,000, three plate lengths 
and varying cylinder to plate spacings. The variations of the 
average Nusselt number as a function of gap spacing was found 
to exhibit an interesting and unexpected trend. The average 
Nusselt number exhibited a minimum value when the cylinder 
was resting on the plate and attained a maximum value at a gap 
spacing of approximately one to two cylinder radii. Flow 
Visualization studies together with measurements of the vortex 
shedding frequency in the cylinder wake indicate that when the 
plate is within approximately one cylinder radius of the 
cylinder, it inhibits the formation of the regular vortex 
pattern that forms in the wake of a free cylinder. This has 
a strong influence on the base Nusselt number datiolala, aay 7RY aeial, 


exerts a controlling influence on the magnitude of the average 


Maeselt number. 
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I. INTRODUCTION 


An experimental investigation was carried out to examine 
ies intluence of wall proximity effects on the heat transfer 
characteristics of a uniformly heated cylinder placed in a cross- 
flow of air. Recently, McComas / 1_/ initiated an investigation 
of the heat transfer characteristics of a uniformly heated 
cylinder attached to a plate and Gnerlich / 2_/ continued the 
same study varying the cylinder to plate spacings. The present 
work is an extensive follow-on study of the work of McComas 
aaa Gnerlich differing in that a range of Reynolds numbers has 
been investigated as well as three different plate lengths. 

The present experimental work was conducted in the Mechanical 
Engineering Department wind tunnel (described in detail in 
Seeger 1fL) in Building 500 at the Naval Postgraduate School. 

Local Nusselt numbers were obtained around the circumference 
of the cylinder at Reynolds numbers of 16,000, 90,000 and 
153,000 without a plate as well as with a plate. In the 
case with a plate, the cylinder to plate spacings were varied 
from d/r = 0 to d/r = 5.33 with L/D = 4 where d = gap width, 
feemeviinder radius, D = cylinder diameter and L = length of 
plate from the centerline of the cylinder to the leading edge 
of the plate (Figure 1). Experiments were also conducted at 
a Reynolds number of 153,000 with L/D = 2 and 8 in order to 
Wee the effect of the length of the plate’s leading edge on 


fee heat transier results. 
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Before proceeding further, a brief qualitative description 
of flow past a cylinder will help to explain and interpret the 
results obtained from this experimental study. Since the 
present experiments examined Reynolds numbers ranging from 
16,000 to 153,000, only the subcritical and critical flow 
regimes will be described here. See Achenbach /3 7 for a 
complete discussion of the various flow regimes. 

In subcritical flow, a laminar hydrodynamic boundary 
layer grows from a minimum thickness at the forward stagnation 
point on the cylinder to a maximum thickness at an angular 
location of approximately 80-85 degrees. At this location, 
fieekinetic energy of the moving fluid is no longer sufficient 
to overcome the adverse pressure gradient present on the 
eyemnder’s surface and the laminar boundary layer separates. 
Upon separation a turbulent wake then develops on the rear 
eeae of the cylinder. 

iimcritieal tioyw, the laminar boundary layer again grows 
from the forward stagnation point to the laminar separation 
memm> which occurs at approximately 80-85 degrees. In contrast 
meme SUbCritical flow case, a transition occurs from laminar 
to turbulent flow. A turbulent boundary layer then reattaches 
megert to the cylinder. This reattached turbulent layer 
eventually separates at a position downstream of the laminar 
separation point. The region between the laminar separation 
Pema and the point of reattachment of the turbulent boundary 


layer is known as a "Separation bubble." 
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The results of the experiments with the cylinder attached 
moeune flat plate indicate a critical flow pattern exists 
mune top half of the cylinder for all Reynolds numbers 
investigated. This flow pattern changed from critical to 
pmeeritical over the top as the cylinder to plate spacing 
Measwanereased to approximately d/r = 0.1 to 0.2. The flow 
over the bottom half of the cylinder revealed critical flow 
immediately after a space existed. Flow over the bottom 
femed trom critical to subcritical in the interval of 
ome d/r = 1.2. 

Plots of the average Nusselt number versus gap spacing 
have shown that the minimum average Nusselt number occurs when 
the plate is attached to the cylinder and the maximum average 
Nusselt number occurs at d/r = 2.0. 

Qualitative studies /4_/, using a flow visualization 
technique in the water tunnel together with measurements of 
the vortex shedding frequency in the wind tunnel, have contri- 
buted a plausible explanation about the trends of average 
Nusselt number curves. 

It is theorized that trends in the average Nusselt 
number curves are a direct consequence of the strength and 


frequency of the vortex pattern which forms in the wake of 


the cylinder. 
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IIT. BACKGROUND 


The present investigation of the heat transfer character- 
istics of a uniformly heated cylinder placed in close proximity 
to a plane surface in an air stream was conducted concurrently 
with an investigation by Géktun / 4 7 of the flow character- 
istics of a cylinder placed near a plane surface. Information 
obtained from the survey of the heat transfer experiments 
Paoved usetul in explaining the nature of the flow field 
observed by Goktun. In the same manner, Goktun's observations 
Made significant contributions to a better understanding of 
the heat transfer characteristics of the cylinder. 

Although extensive investigations have been directed at 
obtaining a better understanding of the heat transfer and flow 
characteristics of “free" cylinders (for extensive references 
on this subject see Cooper, Field and Meyer / 5_/), data 
Pelating to the heat transfer and flow characteristics of 
the cylinder-plate geometry under investigation were not 
egeiiiable in the literature. Such information would be of 
obvious value in the design of heated piping systems placed 
near walls and in the design of heat exchangers. 

As previously mentioned, the present investigation is 
a follow-on study of an investigation initiated by McComas 
[1.7 ané continued by Gnerlich / 2_/. Although McComas 


experimentally determined local and average heat transfer 


Ke 





coefficients around a uniformly heated cylinder attached to 
mee) and Gnerlich did Jikewise for various cylinder to 

plate spacings, both of these investigators found no references 
we theoretical work on this subject. 

Convective heat transfer is significantly affected by 
fitemrttuid flow characteristics occurring adjacent to a solid 
boundary. Since forced-convection heat transfer is a result 
of the movement of fluids and the mixing of fluid particles, 
the mechanism of fluid flow must be understood in order to 
understand the mechanism of heat transfer and to explain the 
phenomena occurring during the heat transfer process. Because 
Gements, this project was conducted concurrently with Goktun's 
[4/ study of the flow field near the cylinder-plate geometry. 
In explaining the heat transfer results obtained in this 


work, fregquent reference will be made to Goktun’s Concurrent 


flow investigation. 
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Ill. EXPERIMENTAL APPARATUS 


1. Wind Tunnel 

An open circuit, subsonic wind tunnel was used for the 
experiments. The tunnel was specially designed by Kenney 
Engineering Corp., of Monrovia, California for research 
quality work. The entire exterior dimensions are given in 
Figure 2. The tunnel has a test section which is 20 inches 
wide, 28 inches high at the entrance and expands to 31 inches 
Euene EX1t bo provide for boundary layer growth. It is 
powered by a 75 HP, 1750 RPM variable speed controlled D.C. 
Moecor Which provides a belt drive to an adjustable pitch 
mameixtial fan. The fan is 45 inches in diameter and is 
eaoagoke of Supplying air at 10,000 to 70,000 cfm, producing 
Wind velocities ranging from 10 mph (14.7 ft/sec) to 175 mph 
(256.7 ft/sec). The tunnel is equipped with six graded mesh 
screens to control turbulence. A 24 inch water manometer 
connected in to the stilling chamber and test section may be 
used to measure velocity. Utilizing a hot wire anemometer, 
a turbulence intensity level of 0.2% was found in the wind 
firme! at a free stream velocity of 98.1 ft/sec. A mercury 
thermometer was installed in the stilling chamber to measure 


the free stream temperature. 
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wae Liguid Crystals 

Microencapsulated cholesteric liquid crystals were used to 
obtain both qualitative and quantitative heat transfer and 
fluid flow information on the uniformly heated cylinder. The 
miquid crystals provide a visual display of temperature over 
discrete, reproducible temperature ranges (called the "event" 
temperature range) by means of a change of color (see Cooper, 
Field and Meyer / 5_/). The visible color Spectrum, which 
falls within a certain temperature band, has an accuracy of 
oer, All experiments in this investigation were conducted 
Using a liquid crystal designated as NCR S-43 having an event 
temperature range of 

Red : 109.8°F (43.2°C) 

Green: 111-0°F (43.9°C) 

Bie: 112.3-F (44.6 -C) 


3. Temsheet 

An electrically resistive carbon impregnated paper known 
as “Temsheet" was used as the heating element in these studies. 
Temsheet is a thin, highly flexible, fibrous material contain- 
ing no wires or ribbons and therefore is free of localized hot 
peos, burnout, and breakage problems. The nominal thickness 
foo.) cm and the electrical resistivity is about 2.5 ohm-cm. 
Using the Temsheet it is possible to heat large areas within 


fees) than 2% deviation from point to point. 


vad) 





4, Cylinder 

The test cylinder utilized in obtaining the heat transfer 
data was the one used by Gnerlich / 2.7. It was made of three 
cylinders attached end to end (see Figure 3) having a total 
length of 20 inches, corresponding to the test section width, 
and a diameter of 3 inches. The two end sections were hollow 
wood cylinders connected by a wooden brace, curved to match 
pmescylinder contour. Qheweemmral cylinder, henceforth 
designated as the test section, was formed from a sheet of 
Temsheet. It was tightly packed with glass wool which served 
the dual purpose of providing insulation against possible 
internal heat losses due to conduction and free convection 
and also provided the paper cylinder with added resistance 
to deformation. The circumference of the Temsheet cylinder 
Was marked at five degree angular intervals in two locations 
Peacies apart for accWeavemmeeation of the liquid crystal 
isotherms during the experiment. Aluminum foil conducting 
tape served as electrodes. The aluminum tape was circumfer- 
entially attached on the inner surface of the Temsheet. Siiver 
conducting paint was applied to the interface between the tape 
and Temsheet to insure good electrical contact was maintained 
between the two surfaces. Two lead wires were connected to 
m@emaluminum foil througmeememend of the cylinder. The lead 
wires were attached to a regulated DC power supply (model: 
LAMBDA LK 3454 FM). The resistance of the Temsheet was 


measured by a digital multimeter (model: FLUKE 8100A). 
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fa Frlates 

The plates used in the experiments were constructed of 
plexiglass and were 0.5 inch thick, 20 inches wide and had 
total lengths of 40, 46 and 58 inches, respectively. The 
plates were installed horizontally in the wind tunnel in such 
a way that the leading edge of the plates extended 2, 4 and 
8 cylinder diameters (i.e. 6, 12 and 24 inches, respectively) 
ahead of the cylinder centerline. The plates were beveled 
on the under side of the leading edge with an angle of 12 
degrees. The plates were capable of being located at any 
position relative to the cylinder (Figure 4), from attached to 
a maximum distance of 8 inches away (i.e. d/r = 5.33). The 
plates were moved vertically by means of a threaded rod 
extending through the tunnel floor, and a rear brace which 


locked in position by means of a split nut. 
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IV. EXPERIMENTAL PROCEDURE 


The present investigation was conducted using the same 
test facilities as employed by Gnerlich / 2.7. Gnerlich, 
however, performed experiments only for a Reynolds number of 
153,000 (based on the cylinder diameter) and a plate leading 
edge of 12 inches. In the present work, experiments were 
Conducted at Reynolds numbers of 16,000, 90,000 as well as 
153,000 and the plate leading edge was varied from 6 to 24 
inches. 

To initiate the experiment, a sufficient voltage (from 
iZeonl5 volts) was applied to the cylinder to insure that the 
liquid crystals reached their event temperature and changed 
color over the entire surface of the cylinder. This was done 
to check for hot and cold spots which would indicate Aone 
packing of the glass wool inside the cylinder. The wind 
tunnel was then started and set to a predetermined velocity 
at which it was maintained throughout the run. The applied 
voltage was adjusted until the first faint red color trace 
was observed. The power supply was then turned off and the 
multimeter was switched from DC voltmeter to ohmmeter. The 
same procedure was also followed at the end of the experiment. 
imernias manner, an accurate value of the paper resistance, 
which was taken to be the arithmetic mean of initial and 
final resistance values, was determined. The multimeter was 


then switched back to the DC voltmeter and power was re-applied 
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to the test section. A steady state condition was usually 
attained in from five to ten minutes depending upon the 
wind velocity. Generally, at least 40 minutes were spent 
to take the first data point in order to be extremely certain 
about steady state conditions. 

ihe first color band observed om the cylinder was at the 
location of the value of the lowest local Nusselt number. 
mites voOltagse was then increased such that the liquid crystal 
color bands would appear at various angular positions around 
the circumference of the cylinder. The temperature obtained 
in this manner could be compared to the air stream temperature 
which was obtained from a thermometer mounted in the stilling 
chamber. Using this temperature difference and power dissi- 
pation, the local Nusselt number based on the cylinder diameter 
could be obtained around the cylinder (see Appendix A). 

Due to the method of construction, there was a seam along 
the length of the test cylinder. Measurements within 30 
degrees of a sector of the seam were considered unreliable. 
The majority of Nusselt number data were collected with the 
seam located at 180 degrees. The Nusselt number data in the 
vicinity of 180 degrees were taken by turning the cylinder 
Suen that the seam was facing forward at 0 degrees. several 
other seam positions were also checked but it was ultimately 
@ecrded that the 0 and 180 degrees locations were the best 
for the seam position. In all experiments, the location of 


emo ligutd crystal color transition points were easily 
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distinguishable around the circumference of cylinder except 
in the region of 0 to 45 degrees and 160 to 200 degrees. 
Especially in the forward stagnation region the temperature 
gradient was very shallow and the whole region appeared 
memnermal. in the vicinity of the rear stagnation point, 
the temperature gradient was indistinct due to the turbulent 


behavior of the flow in the cylinder wake. 
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Y. RESULTS 


The results shown in this section were obtained from 
the heat transfer experiments discussed in Section y~T. The 
Beeults are classified in two groups: 
A. LOCAL NUSSELT NUMBERS, Nu 

a. Typical data and a sample calculation for the heat 
transfer coefficient and the local Nusselt number are shown 
in Appendix A. 

fee teures 5 through 12 show the variation of the local 
Nusselt number with gap size at a Reynolds number of 153,000 
and a plate leading length of 12 inches (L/D = 4). These 
plots also show the variation of the pressure coefficient 
as a function of gap size. All pressure data were taken 
from Goktun's / 4 / results and are informative in that they 
show the relationship between surface pressure and heat 
transfer characteristics. 

ec. Additional plots of the Nusselt number versus gap 
size for Reynolds numbers of 16,000, 90,000 (with L/D = 4 ) 
and 153,000 (with L/D = 2 and 8 ) can be found in Appendix D, 
Figures 24 through 55. Local Nusselt number distributions on 
the surface of a free cylinder are also contained in Appendix D, 


Heemres 56 through 58. 


B. AVERAGE NUSSELT NUMBERS, Nu 
a. A graph of Nu versus d/r (the cylinder to plate 


spacings) for Reynolds numbers of 16,000, 90,000 and 153,000 
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and L/D = 4.0 is shown in Figure 13. A comparison is made 
with Gnerlich's / 2_/ results on the same graph. 

b. A graph of Nu versus d/r for Reynolds number of 
lMee,000 and L/D = 2, ¥ and 8 is shown in Figure 14. 

c. A graph of normalized drag coefficient (Ca/Cy_ ds 
Strouhal number (S/S. ) and average Nusselt number (Nu/Nu,, ) 
versus d/r is shown in Figure 15. The “co" subscript refers 
fo values obtained with no plate in the tunnel. 

d. A graph of normalized Strouhal number and local 
base Nusselt number eee ieee) versus d/r is shown 


in Figure 16. 
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Figure 5. Local Nusselt number and pressure coefficient on the 


surface of a cylinder placed near a plane surface for 


bBew= 153,000, L/D = 4, d/r = 0.0. 
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Figure 8. Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 


hom ence, L/De= 4, d/re= 1.0. 
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Figure 9. Local Nusselt number and pressure coefficient on the 


surface of a cylinder placed near a plane surface for 


Ron wlenO00 Daw, d/r = 2:0. 
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Figure 10. Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface for 


Rees l>oneO0 st) = Yaad /r = 3.0. 
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Figure 11. Local Nusselt number and pressure coefficient on the 
surface of a cylinder placed near a plane surface 
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Figure 12. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Siimiaccmncrem helm-9155,000, L/D = 4, d/r= 5.33. 
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Voi CLS CuSO LON 


Because of the relationship between the flow characteristics 
(i.e. pressure distribution around the cylinder, pressure drag 
and vortex dynamics) and the heat transfer characteristics of 
the cylinder, the parallel flow study of Goktun / 4 7 will 
be referred to in some parts of the following discussion. 
Goktun obtained his results using an unheated aluminum cylinder 
placed in proximity to the plane surface. 

me feneral nature and trends of the results obtained 
from the present investigation can be discussed in two 


groups: 


Pee LRENDS IN LOCAL NUSSELT NUMBERS, Nu 

At d/r = 0 

The initial heat transfer measurements were made with the 
Maatormly heated cylinder attached to the plate at a Reynolds 
member of 153,000 with L/D = 4. The data obtained in this 
arrangement agreed within the estimated experimental uncertainty 
of 5% (See Appendix B) with the results of McComas / 1_/ 
and Gnerlich / 2_/. A local maximum and minumum value of 
the Nusselt number were once again observed near the point 
of attachment of the plate and cylinder and are indicative of 
a large trapped vortex as proposed by McComas We da Jee in the 
flow visualization experiment, at d/r = 0, Gdktun / 4 / also 
observed that a series of small, circular, three dimensional 


vortices were also present in the upstream cavity as well as 
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a large turbulent vortex. The large turbulent vortex occurs 
within the separated region where the flow separates from the 
ieu plate and reattaches to the cylinder. The point of flow 
reattachment occurs at approximately 140 to 170 degrees from 
mee point of attachment of the plate and cylinder. ‘he 
existence of this large turbulent vortex produces the observed 
local maximum and minimum values in the Nusselt number curve. 

The existence of a separation bubble on the top of the 
cylinder was the obvious mark of critical flow. A laminar 
boundary layer started to grow on the forward half of the 
cylinder and separated between 80-85 degrees and a subseguent 
ieaiettion to turbulent flow occurred resulting in the 
reattachment of a turbulent boundary layer at about OO Eales” . 
The reattached turbulent layer ultimately separated at 
approximately 130-135 degrees. 

The decrease in the local Nusselt number (see Figure 5) 
at the aft end of the cylinder was unexpected. The flow 
Visualization studies clearly showed, however, that when the 
cylinder is resting on the plate a large wake with very little 
fluid motion develops on the downstream side of the cylinder. 
The absence of fluid motion therefore explains the decrease 
famusselt number on the back of the cylinder. 

feed / i = 0.25. tod/r = 95.353 

Emether experimental results at a Reynolds number of 
153,000 were obtained with the plate as far removed from 


the cylinder as possible, d/r = 5.33. 


8 


Examination of the polar plots in Figure 6 to 12 for 
Significant trends or unusual nature show two items of 
interest: 

1. A shift in the location of the forward stagnation 
point, rear stagnation point and the laminar separation [evo icles 
(see Figure 17). 

Starting from d/r = 0.25 and moving the plate to 
d/r = 2.0, it appears as though an imaginary axis initially 
passing through 340 and 160 degrees is pivoted about the 
origin. As this axis turns clockwise through 20 degrees of 
arc to 360 degrees, the angular position of the stagnation 
and separation points changed correspondingly up to d/r = 2, 
and then reached a constant behaviour (see Figure 6 to 9). 

iigosmobserved shitting of slagpnation polmss wac 
definitely not affected by the variation of the length of 
the plate's leading edge and/or the Reynolds number as shown 
in Appendix D, in Figures 40 to 55. 

Waele transition from critical to subcritical flow. 

ie Teatiet tron heOmscritlocal GO suber. gical imilow wae 
investigated in order to find the d/r ratio associated with 
fmmenpransition on the top and bottom of the cylinder. At 
a Reynold's number of 153,000 and L/D = 4, the flow was 
critical over the top of the cylinder with the plate attached 
pag was subcritical over the top at d/r = 0.25. By alternately 
fearying the plate to cylinder spacing, carefully adjusting the 
molLcvace and noting the liguid crystal display, the transition 
Ppolnt from critical to subcritical flow was precisely found 


meed/r = 0.1. 
Why 
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A transition from critical to subcritical flow over 
the bottom of the cylinder was found at approx imere Ly 
d/r = 1.2. In this case, although the pattern of the 
separation bubble displayed by the liquid crystal band was 
less pronounced than over the top, it was clearly observed 
momear as d/r = 1.2. 

As also noted by Gnerlich / 2_/, flow over both the 
top and bottom of the cylinder are definitely critical for 
the stated d/r ranges but there is a significant difference 
mieune degree of criticality on the top and bottom of the 
cylinder as far as the magnitudes of the local Nusselt numbers 
pee concerned. This is seen in Figures 5 through &. 

Further study was also conducted in order to see the 
influence of the length of the plate's leading edge and 
Reynolds numbers. The results are shown in Table 1. Table 1 
aeons wine transition points from critical to subcritical flow 
that occurred over the top and bottom of the cylinder as a 
function of Reynolds number and L/D ratios. It is worthwhile 
femeote all points fell within the range of 0.1 = d/r ©£0.2 
meteethe top. The transition points for the bottom were 
fepmd in the interval of 0.8 © d/r = 1.2. Also to be noted 
imebable 1 is that the transition from critical to subcritical 
flow at various d/r ratios is apparently independent of 
L/D ratios but is mildly dependent on Reynolds number. As 
the Reynolds number increased, the d/r ratios for transition 
on the top of cylinder decrease slightly; on the bottom of 


the cylinder the behaviour was reversed. 
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mee 1, Transition from critical to subcritical flow as a 


function of Reynolds number, d/r and L/D ratios. 


B. TRENDS IN AVERAGE NUSSELT NUMBER, Nu 
In order to gain further information about the cylinder 
Meee transter characteristics, the average Nusselt number, 


Nu was computed from the formula 


Za 


ee 
Nu = aa f Nu:de 


O 


The integration was carried out on an IBM-360 computer using 
Simpson's Rule of Quadrature with a step size of 5 degrees. 
Figure 13 shows the variation of the average Nusselt 
number as a function of gap spacing for three Reynolds 
number and L/D = 4. 

A comparison of the results obtained at ne — lbs ,O0er and 
L/D = 4, with those obtained by Gnerlich indicates that the 
trends are essentially the same. The sero rig. Neat tee 


d/r = 0.5 in Gnerlich's results was not observed in the 


aa 





present experiments but, considering the magnitude of the 
decrease and the 5% experimental uncertainty involved in 
both experiments, it is not felt to be Significant. The 
maximum Nu was observed at d/r = 2.0 and it is readily 
Emeparent that the presence of the plate radically altered 
the magnitude of Nu in the range 0 © d/r © 2. Beyond a 
d/r of 2.0, Nu approached a nearly constant value. 

it is important to emphasize that the trends of the Nu 
versus d/r curves are nearly the same for all Reynolds 
numbers. 

In order to see the effect of the length of the plate's 
leading edge, results were obtained for a Reynolds number of 
153,000 for L/D ratios of 2, 4 and 8. The results, which are 
shown in Figure 14, clearly indicated that there is no 
significant effect on the trends in the Nu curve due to 
changing the L/D ratios. Godktun / 4 / also noted that 
changing the L/D ratio had no significant influence on the 
teayszeand lift characteristics of the cylinder. 

For the sake of simplicity, the flow characteristics 
(pressure drag and vortex shedding frequency) and the heat 
transfer characteristics (average and base local Nusselt 
number) are combined and plotted in normalized form in 
Figures 15 and 16. 

Goktun's qualitative flow visualization study / 4 / in 
the water tunnel together with the vortex shedding frequency 


measurements in the wind tunnel indicated that vortex dynamics 


control the drag on the cylinder. Figure 15 clearly indicates 
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the interrelation between the drag and heat transfer. It is 
seen that as the drag coefficient increases, the average heat 
transfer coefficient also increases. It is therefore concluded 
that the vortex dynamics are responsible for the trends 
observed in both the drag and average Nusselt number curves. 

In the range 0 Cd/r(1, Gdktun found that the natural or 
regular vortex pattern formed in the wake of a free cylinder 
is inhibited by the presence of the plate. The plate Slgenifi- 
cantly reduces the vortex shedding frequency. This in turn has 
an effect on the base pressure (see Figures 5 through 8) and 
base Nusselt number (Figure 16). The more stagnant the wake, 
the higher the base pressure and the lower the base Nusselt 
number. Therefore, the stagnant wake results in poor heat 
transfer. This trend, in turn, influences the drag and 
average Nusselt number (Figure 15). 

Goktun / 4 _/ observed that strong and violent vortices 
were formed behind the cylinder when d/r was approximately 
1.0 to 2.0. These vortices transfer high-momentum cool fluid 
from outside the boundary layer in toward the rear surface 
of the cylinder. The cool fluid in the vortex which impacts 
on the hot cylinder surface has a large capacity to remove 
energy. The "scrubbing" motion of these vortices building 
up behind the cylinder and breaking away to be carried down- 
stream in the wake, results in high rates of heat transfer 
and this results in a large base Nusselt number. This trend, 


iieeturn, influences both the drag coefficient and the average 


Nusselt number. 


LO 





pestne plave is meved iurther from the eylinder, the 
drag and Nusselt number approach their "free" cylinder 


Values. 


59 





~ Va a | Oe 


" ee 
a a. a - 
abr i '40 Fame ad ae es . 






Voie ow ary 


This experimental investigation was carried out to gain 
an understanding of the heat transfer behavior of a btalsLameye aly 
heated right circular cylinder placed near a plane surface. 

A concurrent study by Goktun / 4 / was also carried out to 
eumey the lift and drag characteristics of a cylinder in the 
same configuration. Local Nusselt numbers were measured around 
imesecylinder circumference for various plate to cylinder 
spacings, plate lengths and Reynolds numbers. The average 
Nusselt numbers were obtained from the local Nusselt number 
meets butions by numerical integration. 

An unusual trend which was observed in the plots of both 
the average Nusselt number and overall drag versus cylinder to 
plate separation curves lead to a flow visualization experiment 
and hotwire anemometer measurements of the vortex shedding 
frequency in the cylinder wake. Based on these observations 
and measurements, it was concluded that the structure of 
vortex pattern that develops in the wake of the cylinder 
exerts a controlling influence on the shape of the average 
Nusselt number and drag curves. In the range of a plate to 
cylinder spacing of one cylinder radius, it was found that 
the presence of the plate inhibits the formation cf the 
@epular vortex pattern that forms in the wake of a free 


cylinder. This significantly affected both the cylinder 


oes 





base pressure and local Nusselt number distribution which, 
in turn, influenced the drag and average Nusselt number 


eurves. 
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APPENDIX A 
HEAT TRANSFER DATA REDUCTION 


The procedure followed in reducing the heat transfer 
eeiea 1c Shown below. in this manner the local Nusselt 
moe rs at various locations were obtained for the case 


Sree uniformly heated (constant heat flux) cylinder. 


HEAT FLUX 


The heat flux produced by Joulean heating is given by; 


H 


where: E Voltage applied to the test section, Vols), 





Ro = Ova leeecicatance Oi Lemsneect, / ohms_/ 
A = Total heated area of Temsheet ( =0.315) [tt*7 
2 
q, = Rate of total heat transfer, V 
Q- £t* 
| Ve BIU 
Mem@yersion: 1 ———5 = 3.413 ne oce 
Osert heait 
a °)|CWW 00925 R R heo-ft 
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HEAT TRANSFER COEFFICIENT 


The total heat transfer coefficient is determined by 
dividing the heat flux by the measured local temperature 


Mieererence (T - TL, ): 


wae qd a Gy 
a a ee rN 
ni = he + Phe 

qwy 

A, AD ny. 

= nate © 2 
Dy i i» I) eae 

where: 

= _ Ppt Ew 

ny =— Tet Veetleat teanorer COeCtTiLICclent, a a 
hr-ft°- ~F 

T = cylinder surface temperature at angular 


location @, as determined by liquid crystals, 


To = air stream temperature, F 
: ee fs Ba 
Ie = convective heat transfer coefficient, 
: he = aye 
: : Pe rev 
Wa radiation Neat wransiter COertili cient, 5 
1% 
hr - ft - 
_ Bak 
¢ = Stefan-Boltzman constant = 0.1714 x 10 oy 
hies= CG 


G12 radiation exchange factor between cylinder 
and surroundings = Oss, 
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NUSSELT NUMBER 


cylinder diameter = ia ri 


where: D 


Kp =~ cleuialwcOldUe miivemny Olea ber wai ane amet 


[Sete ee 
film temperature Tp =e —_-2E 


NUMERICAL EXAMPLE: 
A sample calculation of the local Nusselt number is 


shown below using the data (see page 57 ) obtained at a 


2 hr - ft - CF 


fewmields number of 153,000, @/r = 4.0, L/D = 8 and angular 


location of @= 63 degrees. 
V = 31 volts 
Pee ks 3 
T = 111°F = 571 °R 
T= 70°F = 530°R 


MW=T-T = 41°R 


feo ~ 9-9 


meee x 108 ee 





- i 
D= 3 = oR 
Z 
Qw as 160) ee R ay 10.842 eas 7 Sos, — ree 
no = OH _ 783.39 2 ag. BW 
G AT Z O 
a hie it oS 





Nu 


Nu 


—— 


I 


i 


iI 


Aen Dam y(n a TS ) 


(9.1714 x 107°) (0.9) (572 + 530) (571% + 530%) 





1 leu) 
Zz O 
Nie iG =) FF 
h, -h,= 19.1 - 1= 18,1 ——#2— 
eee tt re 
legs velo 
Ie O 
i SOc. 
e 
Bee : 
0.0154 A (evaluated at film temperature 
hee ete ae Ono, SE) 
1 
Oe Selon) Gee) 
Ke 0.0154 


294 at 6= 63 degrees. 
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DATA 


[rs 4.0 
L/D= 8,0 


Ra = 153,000 


Ug@ 98.1] ft/sec 


in-il20 


a arte 


Rp = 13.3 Ohn. 
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UNCERTAINTY ANALYSIS 


Since there are uncertainties in the measured quantities, 
an uncertainty analysis was performed to ascertain the degree 
of uncertainty in the experimental results at three ranges of 
Reynolds numbers of 153,000, 90,000 and 16,000. The method 
of Kline and McClintock / 6_/ was used. 

Being conservative, the following uncertainties were 


Setimatea in the measured variables: 


miolttace, V = 0,1 volts (20:1) 
resistance, R = 0.4 ohms (208 T 
temperature, T+ 1.0 °F (20:1) 
area, At 0.001 £t* (20:1) 
mrameter, D = 0.001 ft (207) 
wall manometer, H = 0.05 in - Ho (20:1) 
radiation exchange 

Pac tor de = 0.05 (20:1) 


The uncertainties associated with the properties of air 
were considered to be negligible. 

The uncertainty analysis is performed for the data 
obtained at Reynolds numbers of 153,000, 90,000 and 16,000 
as follows; 


i Lllustration of procedure of uncertainty analysis for 
me = 1559, 0007: 
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Total Heat Transfer Coefficient, hy 3 
As shown previously in Appendix A, the total heat 


transfer coefficient was calculated from the equation: 


ihe uncertainty in the heat transfer coefficient is then 


expressed as: 











For the run at a Reynolds number of 153,000, L/D = 8, 


67/r— 4 and an angular location of 63 degrees: 








V = 31 0.1 volts 

R = 13.3 = 0.4 ohms 

pe = 4 = 2.0 CF ( T = 111°F = 571°R, Tea = 70 F = 530 R) 
fee = 0.314 = 0.0013 £t- 

H E220) = OO Sain 50 
ome = 19.1. a - 

hele en = 
W 
Ono or Qe =) eas BTU 
n ° Zz O 
ne ‘ eS Eto wee 


This result does not take into consideration the 
uncertainty in the angular location of the liquid crystal 


: ; + 
band. This uncertainty is estimated to be - 3 degrees. 


De 





The uncertainty in the total heat transfer coefficient 
remained relatively constant around the circumference of the 
cylinder. However, the uncertainty in the angular location 
of the liquid crystal band varied. An estimated maximum 
uncertainty of = 5 degrees occurred in the forward and rear 
stagnation regions and a minimum of = 1 degree occurred in the 


regions near the separation and reattachment points. 


meacdiation Heat Transfer Coefficient, ee 





The equation used to calculate Mk was: 


2 Z 
pe CAMO pea ean) Tree. |) 
es gil 
defining T, = —- (T+ Te) 
fet 3 
= ~ oO: ° 
ny 4 ee Ph 


eerqe]e Go 1S constant, it is not considered in the uncertainty 


Baetysis. The uncertainty in ee is then expressed as: 











3w 2 
UW) 
me i 
n. om 
U) 
h 
BLU 
ae = 0.06 or w, = 0 06 5 
ig ie hye = tb oe 
i 
h. = 1.00 + 0.06 = 
= Meee Fe 2 Op 
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iocal Nusselt number, Nu: 


lige VhOocalaNusselb number Was Galculsated from the 


equation: 
D °D 
Nu = 
Ke 
Therefore, 
Zz 
SN 
Nu 











Since the uncertainties in the properties of air are considered 
to be negligible and the uncertainty in the cylinder diameter 


Mem@eclicibly small; 





Siu. os a 
———- 0.044 or uw, = lies 
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Reynolds Number, Re : 


The equation used to calculate Reynolds number was: 





UW, oe 
Leo 
U 
inerefore 
a 
Re 2 6; 
Re UO; 





Since the uncertainty in the kinematic viscosity and in 


the cylinder diameter are negligible: 


U) 


MiemevaitGity Uj. is proportional to the square root of 


the wall manometer reading, hence: 


"Re _ 1 _‘H 
- 2 H 
” Re 
= 0.011 or ®,. = 2000 
R = 
S 
Re = 153,000 * 2000 


meine uncertainty analysis for the run at Re = 90,000: 
For the data obtained at a Reynolds number of 90,000, 


imo= +, d/r = 2.0 and @= 67 degrees: 
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Uo = 28.0 2 Gy Welty 

R = 13.0 + 0.4 ohms 

teen oie ( Pri = "57k T= 7.4 or = 
A = 0.314 t 0.001 ft” 

He ="0n76= 0.05 in = H,o 


Z 
eimilarly, the uncertainties calculated in following 


wartaples were: 


U) 











h 
t= 9,04 or u =t0.7 ——BWU 
t + hy - £t© - OF 
fe, =) G6) =o BIU 
t 2 O 
hr - £t° - OF 
Oh 
‘= 0.06 or w = 0.06 = 
e r hr - ft*© - OF 
h, = 1.00 t 0.06 aS 
hr - ft* - CF 
“Nu 
— 0.044 or On = dust 
Nu = 259 * 11 
eRe 
Sn eee ee 28 


hey 790,000 +) 3000 
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532.4°R) 





oe Jhe uncertainty analysis for the run at Re = NG OC Ons: 


For the data obtained at a Reynolds number of 16,000 
L/D = 4, d/r = 2.0 and 9 = 65 degrees: 


V 21.1 £0.1 volts 


R 13.4 £ 0.4 ohms 


Hi 


AT = 43 $1 gi Me IP = ei ws ben a Bose) 


a Z 


OO OO) iat 


A 


The following uncertainties were calculated: 








De i 
eer 0-029 or & == 9.3 oe 
t hr - £4° — °F 
ie = Sos Bales oa 
t a O 
hr - ft° - OP 
hh 
—- = 0.06 orw, == 0.06 my 
r r nr - £t° - OF 
h_= 1.00 © 0.06 BIU 
In 2 O 
ne ah a 1 
GU) 
—NU_ = 9 on or On 2 
1B 
Nu 


For this particular Reynolds number it was impossible to 
measure the free stream velocity using the wall manometer. 
The hot wire anemometer was therefore used to measure the 


velocity. Thus, the uncertainty estimated in U,, was: 


6ly 





lee = 2 ft/sec (20:1) 





Re = 16,000 = 3000 
The high uncertainty in Ue and consequently in the Reynolds 


number was primarily due to the unstable operation of the 


fan motor at its lowest operating point. 
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APPENDIX C 
WIND TUNNEL VELOCITY PROFILES 


For the three ranges of Reynolds numbers of 16,000, 
90,000 and 153,000 and with L/D = 4, velocity profiles were 
obtained in the tunnel test section (see Figure 18 through 
23). Data were collected at d/r = 0 and d/r = 5.33 (see next 
Maes for data). 

The hot wire anemometer was used to probe the test 
Peeion at various positive and negative Y - locations at 
90 degrees and 270 degrees, measured from the forward 
Stagnation point on the cylinder. The velocity profile at 
the test section entrance was also measured at various 
igeatrons and it was found to be fairly uniform. 

The turbulence intensity in the test section was obtained 


for a Reynolds number of 153,000 and was found to be approxi- 


mauely 0.2%, 
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Mreure 18, Velocity profile at the top of the cylinder for a/r 
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Preure 19. 
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Velocity profiles at the top and bottom of the 
eein@er torvd/r = 5.33. 
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| Figure 21. Velocity profiles at the top and bottom of the 
GbintdemmrarwG/re=. 5, 33. 
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Additional plots of the Nusselt number versus gap size 


for Reynolds numbers of 16,000, 90,000 and 153,000. 
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Figure 24. Local Nusselt number on the surface of a cylinder 


placed near a plane surface for 


Re = 16,000, L/D = 4, d/r = 0.0. 


ie 





As) 
£ Af 
< 

AS 
f. 
rN 
(s) 
ven 
© f- 
ee 
Cape 
@ 
Sp 2 - 2 2 2 
fi 6 
abe 
© 
Q\ 2 
ov 
ak Sa 
o ta 
Sak 
“A 
oN 
aX, 
. % 
On , 
ieee 25. 


placed near a plane surface for 


L/D = 4, 


Re 


LOTOOOn 


gO 


Nuss: 


é/r = 


76 


0 BOF 





Local Nusselt number on the surface of a cylinder 








Figure 26. Local Nusselt number on the surface of a cylinder 


placed near a plane surface for 


Roe bonege. 7D = 4) Gr = 0.5. 
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Figure 27. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 


Re = 16,000, L/D= 4, d/r= 1.0. 
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[eure 26. Local Nusselt number on the surface of a cylinder 


placed near a plane surface for 


hee slenOCO, 7D =—+, d/r = 2.0. 
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Figure 29. Local Nusselt number on the surface of a cylinder 
placed near a plane surface for 


en MOOen hy WD = d/r = 3.0. 
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Figure 30. Local Nusselt number on the surface of a cylinder 


placed near a plane surface for 


Rom wllcno oc t/) — 97, d/r = 4.0. 
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Figure 31. Local Nusselt number on the surface of a cylinder 


placed near a plane surface for 


Rememlos000, yea Uaaed/r = 5.33. 
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Figure 32. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Stietace tom ke = 90,000), L/D = 4, d/r = 0.0. 
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Figure 33. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Sine emiommnee= OOMOOO L/D = 4 sd/7r = 0125. 
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Figure 34. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Simmecemnongs Re = 90,000, L/D = 4+, d7r= 0.5. 
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Figure 35. Local Nusselt number and pressure coefficient 


the surface of a cylinder placed near a plane 


surface for Re = = = 


Sener L/D — Yd r = 1.0. 
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Figure 36. Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 


Siemcemior ME Ree=— 908000, L/D = Heed/r = 2.0. 
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Figure 37. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed Nea a lene 


Seco mOn wRoe=90NO0G, L/P = 4, d/r = 3.0. 
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Figure 38. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Sueueecmton Ge = COM000, L/D = 4, d/r = 4.0. 
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Figure 39. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


iieecescmeneCe= OOMO00, b/D = 4, d/r = 5,33. 
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Figure 40. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Seco muome neu loon 000, b/D = 2, d/r = 0.0, 
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Figure 41. 


Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 


VWiecomeonrmnRen=— (5a 0000, L/D = 2, d/r= 0.25. 
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Figure 42. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


SUiGmaccmmonmere =——l5o,000, L/D = 2, d/r = 0.5. 
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Figure 43. 
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Local Nusselt number and pressure coefficient on 


the surface of 4 cylinder placed near a plane 
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Surface for Re = 153,000, L/D eG) i alle Ore 
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Figure 44. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


SUGhICemnomene — 15o,000, L/D = 2, d/r = 2.'). 
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Figure 45. Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 


Siento omuomnee=— Noameode L/D) = 2, d/r = 3.0. 
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Figure 46. Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 


Sree cicemacremmn co — sa ,000, L/D = 2, d/r = 4.0. 
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Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 
surface for Re = 


Slee OOO. L/D = 2,0d/2 = 5)33- 
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Figure 48, 
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Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 
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Figure 49. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Sacco muGmneon= #loemeoO, L/D = 8, d/r = 0.25. 
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Figure 50. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


CUmiccemiomehe = 9 155,000, L/D = 8, d/r = 0.5. 
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Figure 51. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Seeacemigihenee= i5an000, L/D = 8, d/r= 1.0. 
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Figure 52. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


aoc emir emOOOn L/D = 6, d/r = 2.0), 
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Figure 54. 


Local Nusselt number and pressure coefficient on 


the surface of a cylinder placed near a plane 
surface for Re 
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Figure 55. Local Nusselt number and pressure coefficient on 
the surface of a cylinder placed near a plane 


Surtace tor her leemeeO, L/D = 8, d/r = 5.33. 
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Figure 56. Local Nusselt number on the surface of a free 


cylinder for Re = 16,000. 
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Figure 57. Local Nusselt number and pressure coefficient on 


— 


the surface of a free cylinder for Re = 90,000. 
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Figure 58. Local Nusselt number and pressure coefficient 


on the surface of a free cylinder for Re = 153,000. 
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